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Abstract 

Research on interstellar ND3 is reviewed and updated 
from the discovery papers. Results from observations of 
a dozen sources at centimeter and submillimeter wave- 
lengths are presented. The two data sets are consistent, 
but do not constrain the excitation conditions in the ND3- 
emitting gas. The column density ratios of NH3/NH2D, 
NH 2 D/ND 2 H and ND3/ND2H, observed in similar sized 
beams, are «10, and present a challenge for both the 
gas-phase and the grain-surface chemistry scenarios of 
deuterium fractionation. The role of shocks and of CO 
depletion in deuterium chemistry is discussed. 

Key words: ISM: molecules - ISM: abundances 



1. Introduction 

Deuterium-bearing molecules are important as probes of 
the very cold phases of molecular clouds, prior to star for- 
mation. Moreover, the isotopic composition of molecules 
is a valuable clue to their formation mechanism. In the 
gas phase, deuterium fractionation occurs by reactions 
with H 2 D+, CH 2 D+ and C 2 HD+ (Millar, this volume). 
Alternatively, reactions on the surfaces of dust grains may 
enhance molecular D/H ratios. Both methods require low 



the 3.47 /im feature indicates abundances of up to 7% of 
H2O ice (Dartois et al. 2002). However, existence in the 
solid state does not imply formation on grain surfaces, 
since NH3 may just passively accrete after forming in the 
gas phase. Soon after Rodgers & Charnley (2001) pro- 
posed ND3 as test of theories of interstellar deuterium 
chemistry, searches were initiated. 

The rotational energy levels of symmetric top molecules 
like ND3 are labeled by the total angular momentum J 
and its projection on the molecular symmetry axis K. In- 
version splits these levels into states which are symmetric 
(5) and antisymmetric (a) upon reflection in the plane of 
the D atoms. Coupling of the 14 N nuclear spin / with the 
rotational angular momentum J causes additional hyper- 
fine splitting (F = I + J). Current instrumentation can- 
not resolve the even finer splitting due to D, which, ow- 
ing to its small quadrupole moment, is only «200 kHz. 
A complete line list is available from the CDMS catalog 



(Mullcr et al. 2001) at the URL http://www.cdms.de 



Under the conditions favourable for deuterium enhance- 
ment, the strongest lines of ND3 should be the ground 



state rotational line (lg — > 0g) at 309.9 GHz and the 
(1,1) inversion line (1" — > If) at 1589 MHz, which are 
both accessible from the ground. 

Because the Pauli exclusion principle only holds for 
fermions, certain states are forbidden for NH3 but allowed 
in ND3. For example, the ND3 ground state rotational 
temperatures (<20 K) and high densities (> 10 5 cm" 3 ) line also has a lg -> Og component at 306.7 GHz, which 



such as found in dark clouds and pre-stellar cores. The 
grain surface route also needs a mechanism to return the 
molecules to the gas phase where they can be observed 
by the means of submillimeter spectroscopy. The gas- 
phase route can account for the deuteration of HCO + 
and N 2 H+, while grain surface chemistry may be impor- 
tant for H 2 CO and CH 3 OH. 

Although astronomical detections of both singly and 
doubly deuterated NH3 exist (Roueff ct al. 2000), the 
deuteration mechanism for NH3 is not clear. Despite early 
claims, searches for NH3 ice based on the 2.21 and 9.0 fim 



does not exist for NH3. However, spin statistics make this 
line 10 times weaker than the lg — > 0q line, which is why 
it has escaped detection so far. 



2. Submillimeter observations 
Emission in the Jk = 1q — *• On line has been detected 



in the dark cloud Barnard 1 (Lis et al. 2002b), two posi- 
tions in the star- forming region NGC 1333 (van der Tak 
t al. 2002] ), and several positions in the molecular cloud 



features have not been successful ( Taban et al. 2002 ) , but 



L 1689N (Roueff et al, in prep.). Table [l] lists the results of 
a more extensive search, made with the Caltech Submil- 



Proceedings of the conference "Chemistry as a Diagnostic of Star Formation, " University of Waterloo, Canada, 21-23 
August 2002 (C.L. Curry & M. Fich eds.). 
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Table 1. Upper limits for ND3 rotational line emission 



Source 


R.A 


(1950) 


Dec. (1950) 


rms 






h 


m s 







/ // 


mK 


OMC-2 


05 


32 


58.6 


-05 


11 


42 


24 


HH1-C 


05 


33 


51.5 


-06 


47 


57 


25 


NGC 2264 C 


06 


38 


26.3 


09 


32 


18 


39 


NGC 2264 G 


06 


38 


25.7 


09 


58 


54 


31 


Serpens S68N 


18 


27 


15.2 


01 


14 


47 


30 


Serpens FIRS1 


18 


27 


17.3 


01 


13 


16 


28 



lim eter Observa tory. Sources were selected from Shah & 
Woc ttcn (2001) to have high NH2 D column densities, but 



the positions used here are from Jijina et al. (1999) for 
Orion and Monoceros and from McMullin et al. (2000) for 
Serpens. System temperatures were ~500 K and integra- 
tion times 30-50 minutes on-source, using beam switching 
with a throw of 240". For an excitation temperature of 
10 K, the noise levels correspond to column density upper 
limits of 1.4-5.0 10 11 cm" 2 . 



Table 2. Upper limits for ND3 inversion line emission and 
absorption 



Source 


R.A. (1950) 


Dec. (1950) 


rms 




h m s 




mK 


Emission targets 


NGC 1333 


03 26 03.6 


+31 04 42 


11-16 


Barnard 1 


03 30 15.0 


+30 57 31 


9-12 


CB 17 


04 00 35.0 


+56 48 00 


16-22 


L 1400 K 


04 26 51.0 


+54 45 27 


12-23 


L 1544 


05 01 15.0 


+25 07 00 


35-38 


NGC 2264 


06 38 24.9 


+09 32 29 


38-52 


L 134N 


15 51 32.0 


-02 42 13 


15-20 


L 1689N 


16 29 27.6 


-24 22 08 


30-46 


S68 FIR 


18 27 17.5 


+01 13 23 


9-12 


Absorption targets 


W 33 


18 11 19.5 


-17 56 40 


36-46 


W 43 


18 45 00.4 


-01 59 16 


101-120 


37.763-0.216 


18 58 33.1 


+04 07 41 


43-62 


W 49A 


19 07 52.1 


+09 01 08 


70-94 



3. Centimeter observations 



Measurements of the pure inversion spectrum of ND 3 go 
;k to the early days of microwave spectroscopy ( Nuck- 
1 et al. 1953); many lines have been measured since 



(Fusina & Murzin 1994). However, only recently, the hy- 
perfine structure of the (1,1) line near 1589 MHz has been 



resolved ( van Veldhoven et al. 2002 ). 

In 2001-2002, we have used the 100-m Effelsberg tele- 
scope to search for the (J, K) = (1,1) 1589.006 MHz, (2,2) 
1591.695 MHz, (3,3) 1599.645 MHz, and (4,4) 1612.997 
MHz metastable inversion lines. The front end was the 
1.3-1.7 GHz primary focus receiver; the back end was the 
1024-channel autocorrelator. System temperatures were 
18-33 K depending on elevation. Integration times are 
4-10 hours for NGC 1333, CB 17, S68 FIR, Barnard 1, 
L 1400K and L 134N, and 1-2 hours for the other sources. 
Frequency switching with a throw of 1 MHz was used for 
all observations. The telescope has a beam size of 8' at 
this wavelength. Search targets for ND3 emission were se- 
l ected to have jV(NH 2 D) > 3 x 10 13 cm -2 in an 80" beam 
( (3hah fc Woottcn 2001 ). We excluded Orion, where at 
19 cm, the H II region blocks our view of the molecular 
cloud. 

The top part of Table || reports the range of rms noise 
levels obtained for the four lines. Only data unaffected 
by interference with terrestrial signals are listed. For the 
(2,2) and (3,3) lines, «10% of the data were affected by 
interference, but up to 50% for the (4,4) line. In addi- 



tion, OH 1612 MHz emission along the line of sight made 
some of the (4,4) data unusable. Fortunately, the (1,1) 
spectra are clear of interference, and a follow-up search 
with Arecibo may prove successful. 

The inversion lines were also searched for in absorp- 
tion towards bright Galactic H II regions. Targets for this 
search were selected to have A^^CO) >10 1 



measured in 6 cm absorption ( Wadiak et al. 1988 ) and to 
have 20 cm flux densities, taken from the NVSS (45" res- 
olution) of >1 Jy, excluding Sgr B2 because of its low dec- 
lination. The noise levels obtained for these sources (Ta- 
ble |[ bottom part) are higher than for emission sources, 
because the background H II regions contribute signifi- 
cantly to the system temperature. 

The centimeter data limit N(ND 3 ) to < 0.3-1. 5x 10 12 
cm~ 2 , consistent with the submillimeter data. However, 
the inversion lines can not only provide independent es- 
timates of the ND3 column density, they can also con- 
strain its excitation. Column densities derived from the 



Jk 



1 



0q line are uncertain by factors up to 2 



through the excitation temperature. The lower frequency 
of the Jk = If — * If line implies a much lower spon- 
taneous decay rate and virtually ensures thermalization. 
However, being beyond the capabilities of 100-m class 
telescopes such as Effelsberg and the VLA (an effective 
130-m telescope), its detection requires 300-m class tele- 
scopes or larger. 
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proxies of abundance ratios. This approach is valid as 
long as the excitation conditions of all molecules are sim- 
ilar, which is likely true except perhaps for NH3. 

The column density ratios of NH 3 /NH 2 D, NH 2 D/ND 2 H 
and ND3/ND2H are each «10 (Table |). This strong 
fractionation presents a challenge to grain surface chem- 
istry, and requires an atomic D/H ratio in the gas phase 
of ~0.15, which is >10x higher than in the models by 



Roberts & Millar (2000). Unfortunately, observations do 



not constrain the atomic D/H ratio in molecular clouds 
well. While H I 21 cm data indicate H/H 2 ~ 10~ 3 ( [Li k 
Goldsmith 2002j ), observing D I 92 cm is very difficult. 



Radius (AU) 

Figure 1. Model of NGC 1333 IRAS 4 A used to determine 
the ND3 abundance. Based on van der Tak et al. (2002). 



4. Abundance of ND 3 

If the excitation of ND3 cannot be measured, it can be 
calculated through a model of the physical and chemi- 
cal structure of the source. Fig. Ill shows the model that 



To determine the deuterium budget in dense interstellar 
clouds, measurements of the HD 119 /im line, with HIFI 
on Herschel and GREAT on Sofia, will be crucial (see 
Phillips, this volume). 

However, standard models of gas phase chemistry also 
have some difficulty explaining the results of Table 0. To 
reproduce the observed column density ratios, the mod- 
els require deuteron transfer reactions to be much faster 
than proton transfers. Indeed, laboratory data indicate 
that in the dissociative recombination of partially deuter- 



ated ions, H is easier ejected than D (Le Petit & Roucfl 



van I der Tak et al. (20"02l| uscdlo determine the ND3 2002| ). On the other hand, recent experiments on the H 

HD reaction indicate a lower formation rate of H 2 D + 



abundance in NGC 1333 IRAS 4A. The density struc- 
ture is obtained from interferometric dust continuum ob- 



than the value used in most chemical models (Gerlich & 



serv itions, while the dust temperature structure comes Schlcmmcr 2002 ) 



from a self-consistent calculation of the thermal balance 
as a function of radius. These quantities peak toward the 
center of the protostellar envelope. Cosmic-ray ionization 
makes H^, at a rate that does not vary with radius. The 



One potential problem with applying the surface chem- 
istry route to pre-stellar cores is their temperatures which 
are too low for significant evaporation of even the most 
volatile ices to occur. One alter native may be non ther- 
reaction of H^ with HD produces H 2 D + , and since it is mal desorption by cosmic rays (Najita et al. 2001), but 



faster at low temperatures, the H 2 D + distribution peaks quantitative predictions do not exist yet. However, 



at h.rge radii. 



Lis et 



The distribution of ND3 can either be assumed to fol- 
low that of H 2 or that of H 2 D + . These two situations 
may apply to the cases of grain-surface and gas-phase 
formation, respectively. The inferred abundance of ND3 
is 3 x 10~ 12 if ND 3 /H 2 is constant, and 1 x 10~ n if 
ND3/H 2 D + is constant. Maps of the ND3 emission will 
enable us to decide between these two possibilities. The 
Sub Millimeter Array (SMA) will soon provide the nec- 
essary resolution of < 10" at 310 GHz. 

5. Chemical implications 

Table || lists measured column densities of NH3 isotopo- 
mers. Only in the case of NGC 1333 IRAS 4A, abun- 



al. (2002a) propose a connection between shocks and the 



chemistry of deuterium. Indeed, all ND3 detections so far 
have been at or near bipolar outflow sources. Perhaps the 
shocks merely compress surrounding gas, which leads to 
the cold, dense conditions favourable for deuterium chem- 
istry. However, the shocks may also desorb grain mantle 
material, while leaving the grain cores intact. Shock ve- 
locities of ^10 km s _1 would be sufficient. Such shocks 
occur frequently in the interstellar medium. One possi- 
bility is that interstellar gas is periodically shocked by 
ongoing star formation activity. 

Another avenue worth exploring is the connection of 
deuterium chemistry with CO depletion. Freeze-out of 
CO, the major destroyer of H 2 D + , onto dust grains helps 



dan ces of NH3 and ND3 have been derived by radiative 
transfer modeling, so we use column density ratios as 



channeling deuterium into heavy molecules (Roberts & 



Millar 200C ) and correlations between enhanced deutera- 



tion and CO depletion have indeed been observed (Hatchell 
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Table 3. Observed column densities of NH3 isotopomers 



Source 



iV(NH 3 ) Af(NH 2 D) iV(ND 2 H) A(ND 3 





10 14 cm" 2 


10 13 cm" 2 


10 12 cm" 2 


10 11 cm" 2 


K 


Barnard 1 


21 c 


58 6 


48 6 


13 d 


5 


NGC 1333 IRAS 4A 


13. 8 a 


39.0 a 


<4.6 b 


2.9 e 


10 


NGC 1333 D-peak 


19.2° 


37.6 a 


<l.9 b 


5.9 e 


10 


Beam size (") 


37 


22 


22 


25 





a From [Hatchell (2002) 

b Preliminary values based on IRAM 30m observations (Gerin et al., in prep.) 
From Effel sberg observations 

d Data from |Lis et al. (2002b)|; n on-LTE calculation at n(H 2 )=7xl0 4 cm" 3 and T=12 K 
e From 



van der Tak et al. (2002) 



Excitation temperature used to analyze the NH2D and ND2H data 



2002; Bacmann et al. 2002 ). However, in the case of NH3, 
such a trend does not provide direct evidence of forma- 
tion in the gas phase. The formation of NH3 and its iso- 
topomers on grain surfaces requires a high abundance of 
N in the gas phase. While in molecular clouds, most ni- 
trogen is in N 2 flWomack et al. 1992| ), depletion of CO 
suppresses the formation of N2 and enhances N instead 
flCharnley fc Rodgers 2002 ). One key question is therefore 
the major nitrogen carrier in pre-stellar cores. Although 
neither N nor N2 cannot be observed directly, N2H+ can 
serve as a tracer of N 2 . Correlation between the frac- 
tionations of N2H + and NH3 would thus constitute an 
important test of the formation mechanism of interstellar 
NH 3 and ND 3 . 
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